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^ya^ 1 We derive a simple algebraic criterion to select the optimal detector network for a coher- 

^ ' ent wide parameter-space (all-sky) search for continuous gravitational waves. Optimality 

^ ^ ' in this context is defined as providing the highest (average) sensitivity per computing 

^ cost. This criterion is a direct consequence of the properties of the multi-detector T- 

. statistic metric, which has been derived recently. Interestingly, the choice of the optimal 

' network only depends on the noise-levels and duty-cycles of the respective detectors, and 

' not on the available computing power. 

CN ■ 

^ • 1. Multi-detector matched filtering 

' The JF-statistic^ is a coherent matched-filtering detection statistic for continuous 
gravitational waves (GWs). We follow the expressions and notation of our previous 

\ work"'^ (Paper I) on the multi-detector JF-statistic metric. We consider a set of N de- 

' tectors with (uncorrelated) noise power-spectra S'x, where X is the detector index, 

. X = 1, ...,N. Let T be the total observation time spanned by the data to be ana- 

I lyzed. The corresponding multi-detector scalar product for narrow-band continuous 

' waves can be written as 

7: {x\y)=TS-'{xy},, (1) 

5-H ■ 

bX). where boldface notation denotes multi-detector vectors, i.e. {x{t)} = x-^{t). We 

I can allow for the fact that each detector will be in lock only for a duration Tx < T, 

V~j ■ so each detector can be characterized by a "duty cycle", dx = Tx/T < 1. This 

, is a slight, but straightforward generalization with respect to Paper I, and the 

' corresponding noise-weighted time average is defined as 

X 

where the weights wx and the total inverse noise-power are defined as 



00 



iQ)s^ifl^^^ I Q^it)dt, (2) 



5- 



wx = dx^, where S-' = J2dxS^'. (3) 

x=i 

The importance of Eq. ([1]) is that it separates out the scaling with the total obser- 
vation time T and the set of detectors (via S~^), from the averaged contribution 
{xy)g, which does not scale with T or the number of detectors. In terms of this 
scalar product ([T]), the optimal signal-to-noise ratio (SNR) for a perfectly-matched 
signal s{t) can be obtained as 



piO) = VW) = ^TS^\ {^')s- (4) 
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It is obvious from this expression that the SNR increases when increasing the ob- 
servation time T or the number of detectors N . However, here we are interested 
in the case of wide parameter- space searches, in which the highest achievable SNR 
is computationally limited. We therefore need to find the optimal sensitivity per 
computing cost. 

2. Optimizing sensitivity per computing cost 

For simplicity we only consider the sensitivity to an "average" sky-position, so we 
disregard the dependence of (s^)^ to both the sky-position as well as the relative 
orientation of the different detectors. Both should be small effects on average. In 
addition to Eq. (|4]) for the SNR, the second ingredient for the optimal network is 
the computing cost of a wide-parameter search. For the sake of example we consider 
a search for GWs from unknown isolated neutron stars, with unknown intrinsic GW 
frequency /, sky-position a, 5 and one spindown-parameter /. One can show^ that 
in this case the number of required templates Mp scales (at least) as Mp oc T^, 
which severely limits the computationally affordable observation time Tmax- Most 
importantly, however, the number of templates does not scale with the number A'' 
of detectors.^ The corresponding computing cost Cp required to search these Mp 
templates can be estimated as Cp oc NT'' for a "straightforward" computation, 
while it could be reduced down to about Cp (x NT^ if the FFT-algorithm is used.^ 
Generally, we can write 

Cp(xNT'', (5) 

where typically k '--^ 6 — 7 for isolated neutron-star searches. The linear scaling with 
TV comes from the fact that we need to compute the correlation of each template 
with each of the N detector time-series X'^{t). 

The question we are trying to answer is the following: for given computing power 
Cp and a set of N detectors, which (sub)-set of TV < TV detectors {X} C {X} yields 
the highest SNR? Using ([5|), we can express Tmax « (Cp/iV)^^", and inserting this 
into ([3]), we find p(0) cx Cj,^'"^'^^ ■y/7({X}), where the "gain function" 7 is defined as 

N 

7({X})^iV-V'' ^dx5x^ (6) 
x=i 

This simple algebraic function provides the sought-for criterion for the optimal 
detector-network {X}, depending only on the respective noise-floors Sx and duty- 
cycles dx- The optimal detector network is simply the subset {X} of detectors that 
maximizes the gain-function 7({X}). 

This optimal subset can be found in the following simple way: we label the 
detectors X in order of decreasing dxS^^ , and include exactly the first X = 1, N 
detectors in ([6]) where 7 reaches a maximum. It is easy to see that this arrangement 
is optimal, as either adding further detectors, or replacing any term dj^S^^ in the 

sum by another detector X' > iV reduces 7. 
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In the special case of identical detectors, the gain function 7 is strictly monotonic 
with N, and so the optimal network simply consists of using as many detectors as 
possible, reducing the observation time T. 



3. Example application 

As an example, consider a set of "typical" detectors as given in Tabled] The assumed 
parameters are: LIGO (HI, H2, LI) at design sensitivity, with S5 duty-cycles, GEO 
(Gl) at S5 sensitivity, and S4 duty-cycle, Virgo (V2) at design sensitivity, assuming 
a "typical" LIGO duty-cycle. We see that our simple criterion tells us that for a 

Table 1. Example set of detectors with "typical" sensitivities and duty- 
cycles. 





Frequency 


HI 


LI 


H2 


Gl 


V2 






0.71 


0.59 


0.78 


0.97 


0.7 


VSI [10-23/Vh;] 


/ = 200 Hz 


2.9 


2.9 


5.8 


73 


4.4 




/ = 600 Hz 


7.5 
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15 


39 


5.5 
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Fig. 1. SNR gain \/'y{{X}) (assuming k = 6) as a function of the detector network, normalized 
to the single-detector case Left figure: at f = 200 Hz. Right figure: at f = 600 Hz. 



search at / = 200 Hz we should include HI, LI, V2, and H2 for the best all-sky 
sensitivity per computing cost, gaining on average a total of about 40% in SNR over 
HI alone. Similarly, at / = 600 Hz, we find the same set of detectors to be optimal, 
with H2 providing a smaller marginal improvement. 
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